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Abstract 

MHD winds can emanate from both stars and surrounding disks. When the two 
systems are coupled by accretion, it is of interest to know how much wind power is available 
and which (if either) of the two rotators dominates that power. We investigate this in 
the context of multi-polar planetary nebulae (PNe) and proto-planetary nebulae (PPNe), 
for which recent observations have revealed the need for a wind power source in excess 
of that available from radiation driving, and a possible need for magnetic shaping. We 
calculate the MHD wind power from a coupled disk and star, where the former results from 
binary disruption. The resulting wind powers depend only on the accretion rate and stellar 
properties. We find that if the stellar envelope were initially slowly rotating, the disk wind 
would dominate throughout the evolution. If the envelope of the star were rapidly rotating, 
the stellar wind could initially be of comparable power to the disk wind until the stellar 
wind carries away the star's angular momentum. Since an initially rapidly rotating star can 
have its spin and magnetic axes misaligned to the disk, multi-polar outflows can result from 
this disk wind system. For times greater than a spin-down time, the post-AGB stellar wind 
is slaved to the disk for both slow and rapid initial spin cases and the disk wind luminosity 
dominates. We find a reasonably large parameter space where a hybrid star+disk MHD 
driven wind is plausible and where both or either can account for PPNe and PNe powers. 
We also speculate on the morphologies which may emerge from the coupled system. The 
coupled winds might help explain the shapes of a number of remarkable multi-shell or 
multi-polar nebulae. Magnetic activity such as X-ray flares may be associated with the 
both central star and the disk and would be a valuable diagnostic for the dynamical role of 
MHD processes in PNe. 

Subject headings: jets and outflows; accretion, accretion disks; planetary nebulae: 
general; ISM: magnetic fields: MHD; stars: AGB and post-AGB 
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1. Introduction 

Planetary Nebulae (PNe) and proto-Planetary Nebulae (PPNe) are believed to be 
the penultimate evolutionary stage of low and intermediate mass stars (zero age main 
sequence mass range 0.8M Q < M zams < 5M ). PNe appear on the sky as expanding 
ionized plasma clouds surrounding a hot central star. As the resolution of ground based 
telescopes increased, the "typical" shape of a PN went from spherical (Osterbrock 1972) 
to elliptical and/or bipolar (Balick 1987). More recently, deeper and higher resolution 
studies (particularly those with the HST) have shown many PNe with narrow collimated 
features that are better described as jets than bipolar lobes. In some cases the jets appear 
as extremely well collimated bipolar outflows: M2-9 (Balick 2000): Hen-401 (Sahai 2000), 
AFGL 2688 (The Egg) (Sahai et al. 1998). In other cases they appear as FLIERS (Fast 
Low- Ionization Emission Regions), a single knot expanding away from a (usually) elliptical 
nebula. In addition to the presence of jets it also appears that many PNe show evidence for 
multi-shell structures or multi-polarity. In these objects two nested pairs of bipolar lobes 
appear. In multi-shell structures the lobes are aligned along the same axis of symmetry (Hb 
12: Welch et a/.1999) while multi-polar cases show different orientations for the inner and 
outer lobes (M2-46: Manchado et a/.2000). 

While considerable progress has been made in understanding the hydrodynamics or 
magneto- hydrodynamics of shaping bipolar lobes (see Frank 1999 for a review), the origin 
of jets, point symmetry and multi-polar outflows in PNe poses fundamental challenges for 
theory. 

The basic hydrodynamics of bipolar outflows can be described by the classic Generalized 
Interacting Stellar Winds (GISW) model (Kwok et aQ978, Kahn & West 1986, Balick 
1987, Icke 1988). In this scenario a slow (10 km s" 1 ), dense (10~ 4 M yr _1 ) toroidal wind 
expelled during the AGB is followed by a fast (1000 km s^ 1 ), tenuous (10~ 7 M yr" 1 ) 
wind driven off the contracting proto-white dwarf during the PNe phase. The GISW model 
can explain bipolar morphologies and even jets. Point-symmetry and multi-polar bubbles 
however, do not fall easily into line with the concept of a large-scale collimating torus. 

Models invoking a toroidal magnetic field embedded in a normal radiation driven stellar 
wind have recently shown some promise. This so-called Magnetized Wind Bubble (MWB) 
model was first proposed by (Chevalier & Luo 1994) and has been studied numerically by 
(Rozyczka & Franco 1996) and (Garcia-Segura et al. 1999). In these models the field at the 
star is dipolar but assumes a toroidal topology due to rapid stellar rotation. Collimation is 
not activated until the wind passes through the inner shock. Then hoop stresses associated 
with the toroidal field dominate over isotropic gas pressure forces and material is drawn 
towards the axis producing a collimated flow. This mechanism has been shown capable 
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of producing a wide variety of outflow morphologies including well collimated jets. When 
precession of the magnetic axis is included in fully 3-D simulations, the MWB model is 
capable of recovering point-symmetric morphologies as well (Garcia-Segura 1997). It is 
noteworthy that the wind in these models differs substantially from standard MHD driven 
winds in disk and stellar wind contexts, and the fundamental issue of collimation before the 
wind-wind interaction remains to be addressed. In addition it is not clear that such models 
can be applied to multi-polar bubbles due to the difference in collimation orientations. 

Finally neither the GISW or MWB models fully address the origin of the wind. This 
becomes a critical issue when one is dealing with Proto-Planetary Nebulae (PPNe). There 
is growing evidence suggesting that the PN shaping process begins in the PPN phase via 
well collimated jets (Sahai & Trauger 1998). The creation of fast (V > 100 km/s) high 
mass loss rate (M > 10~ 6 M yr _1 ) winds during these early phases when the star is an F 
or G type is not easily explained via classic line driven wind theory as the available power 
is insufficient (Alcolea et al. 2000). Thus for PPN, theorists face a dilemma similar to the 
problem of YSO outflows in that they must explain both wind acceleration and collimation 
(Pudritz 1991, Kdnigl & Ruden 1993, Shu et al. 1994). 

The above considerations motivate further study of magnetized wind models for the 
launching and collimation of outflows from PN. In particular, the basic question of the 
available MHD wind power is important to address. 

Magnetically driven outflows are a leading paradigm for a ubiquity of source classes 
in nature. This includes such contexts as protostellar jets (c.f. Smith 1998; Frank 1999 
for reviews), the Solar wind (c.f. Parker 1979), jets of Galactic accretors (c.f. Mirabel and 
Rodriguez 1999), active galactic nuclei (c.f. Ferrari 1998; Blandford 2000 for reviews), and 
gamma-ray bursts (e.g. Usov 1992; Duncan & Thomspon 1992; Blackman et al. 1996; 
Meszaros & Rees 1997; Ruderman et al. 2000 ). While the parameter regimes are varied 
and the extent and nature of collimation, particle acceleration, and radiation process are 
different in the different settings, an underlying principle remains: magnetic fields can act as 
a drive belt between gravity and energy deposition at large radii, extracting the rotational 
energy of the rotator into an outflow. We note that the potential for accretion disks to exist 
in PPN systems has been raised (Morris 1987; Soker & Livio 1994; Reyes-Ruiz & Lopez 
1999). 

In the case of PN, as in other star-disk or compact object-disk systems, the magnetically 
driven outflows can in principle emanate from both the central object or the disk. Because 
the relative collimation and power evolution of the disk and stellar winds can be different, 
it may not be unreasonable to see outflows of different character coexisting in the same 
object. In order to investigate this properly however, it is necessary to allow that the disk 
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and star are coupled and this requires a model for the disk evolution. 

A careful investigation of the relevant disk formation is presented in Reyes-Ruiz & 
Lopez (1999). The disk in this study forms when a binary system undergoes common 
envelope evolution. The ejection of the envelope involves transfer of angular momentum 
from the secondary to the envelope after which the secondary loses enough angular 
momentum to fall to a separation such that it can fill its Roche lobe and form the disk. 
Reyes-Ruiz & Lopez (1999) discuss a number of important constraints on the constituents 
and properties of binaries which lead to disk formation. The most likely system turns out 
to be an evolved AGB star with mass 2.5 < M a gb/M & < 5 for the primary, with a secondary 
of mass < O.O8M and with an initial binary separation of > 50 — lOOi? (Iben 1991). The 
AGB star will shed 80% of its mass during the common envelope ejection, leaving a post 
AGB stellar core of Mass M c surrounded by a thin residual convective shell. 

The time scale for the disk to move to its inner inner radius is of order or shorter than 
the viscous time scale. Even for systems whose initial outer radii are of order 100-R©, the 
relevant time scale would be no more than a few years. Since the disk would form only after 
the period of common envelope ejection, < lyr, we would expect the disk and the stellar 
shell to form concurrently within a time of order a few years. Because the relevant stellar 
field will be exposed once the envelope ejection occurs, and because any stellar dynamo 
growth time in the disk (discussed later) is less than 1 year, we suggest the rotating stellar 
shell and disk could both support co-operative jets at a time coincident within about a year 
after the disk forms. This represents the initial time appearing in our calculations. 

In section 2, we determine the MHD wind power from PPN stars and disks. This 
requires a calculation of the expected field strength in the disk and a model of how the 
accretion affects the stellar field, the stellar spin, and the disk inner radius. We study the 
cases for which the stellar shell is initially rapidly or slowly rotating, and calculate the 
maximum disk and stellar magnetic wind powers for each case. In section 3 we discuss 
the observational implications of the results of the power calculation and speculate on the 
magnetic shaping. We summarize and conclude in section 4. 

2. Winds from Disks and Stars 

There are many models of magnetic jet production and collimation in the literature. 
The extent and explicit details of collimation are perhaps less agreed upon than the 
principles of launching. Magnetic launching mechanisms can be divided into two basic 
classes, "spring" mechanisms (Uchida & Shibata 1985; Contopoulos 1995 Lynden-Bell 1996) 
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and "fling" mechanisms (Blandford & Payne 1982; Lovelace et al. 1987). In the former 
class of models, the initial driving force is a magnetic spring, i.e. the magnetic field energy 
density is of order the kinetic energy density in the disk and the outflow is driven by toroidal 
field pressure. In the fling mechanism, the initial driving is centrifugal along poloidal field 
lines. Higher up in the corona the wind becomes magnetically driven. In this study we 
compute the total available magnetic luminosity L m from the disk and central star. The 
magnetic luminosity provides an upper limit to the associated wind kinetic luminosity. The 
L m is the integral of the Poynting flux over the surface area of the rotator. We thus have 

L w ~ M w Vl ~ e L m = e J (E x B) ■ dS ~ e J*° \ttR)B p B^RdR, (1) 

where e is an efficiency, and B^ and B p are the toroidal and poloidal field respectively, 
Q and R are the rotational frequency and radius of the magnetized wind source, M m is 
the wind outflow rate, and V w is the wind outflow speed. In what follows we discuss this 
equation in the context of PPN disks and post AGB stars. For all cases of interest, for disks 
the integrand falls off fast enough such that the inner radius matters most. Thus for disks 
we have 

Ldw ~ CdLdm ~ €dB p (R i )B^(Ri)Q d (R i )Ri , (2) 

where fid is the disk angular speed Ld m is the available magnetic luminosity associated with 
the disk, and is the unknown efficiency factor. For stars, we have 

L sw ~ € s L srrl ~ B p {R if )B ( j > {R if )Q if {R lf )R^ (3) 

where £7* is the angular speed of the stellar surface in which the field lines are anchored, 
L sm is the available magnetic luminosity associated with the star, and e s is the stellar wind 
efficiency factor. The above disk and stellar wind luminosities are not independent. The 
poloidal and toroidal fields are coupled through differential rotation in both the disk and 
star, and the angular velocity of the star may depend on the accretion. In addition, the 
inner radius of the disk depends on the stellar magnetic field. Thus (0) and (ff) are very 
much interdependent. 

We now proceed to solve for the coupling relations and consider the stellar and disk 
MHD wind luminosities for two cases: 

• The star is initially a strongly magnetized rotator independent of the accretion. 

• The star is initially a slow rotator but is spun up by the disk accretion. 

As we will see, the evolution in both cases depends on estimates of the disk and star 
magnetic fields. We must also solve for the disk inner radius, and the stellar surface angular 
speed evolution. 
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2.1. Magnetic Field of Disk 

If the jet/wind is due to a magnetic field in the disk, this field must be produced in situ 
for the PN case. This is because the companion star, the progenitor of the disk, would have 
to be a brown dwarf or a Jupiter type planet (Reyes-Ruiz & Lopez 1999) which typically 
would only have fields of order or less than 10G. Even a field several orders of magnitude 
larger would be far too small to have any influence on the disk, which forms from shredding 
this mass, if the disk field were only that resulting from flux freezing. We estimate the 
field strength roughly from dynamo or turbulent amplification in the disk. An alternative 
dynamo approach to ours is given by Reyes- Ruiz & Stepinski (1995). 

For an accretion disk whose angular momentum is transported by a magneto-shearing 
instability (c.f. Balbus-Hawley 1991, 1998), the time scale for growth of unstable modes 
and the time scale for the largest eddy turnover time are both approximately equal to the 
rotation time scale. Very roughly, the eddy turnover time can be written as L/vt where vt 
is the dominant turbulent velocity and L is the dominant correlation scale. We can estimate 
the strength of the mean magnetic field by first estimating the total magnetic energy. Using 
the Shakura-Sunyaev (1973) viscosity prescription 

v = a ss c s H ~ v T /Q, (4) 

where a ss ,c s ,H and L/vt are the disk viscosity parameter, the sound speed, the scale 
height and the dominant eddy turnover time respectively. Using the fact that turbulent 
stretching leads to Va ~ i>t 3 where Va is the Alfven speed (ignoring a possible factor of 
on vt, e.g. Balbus & Hawley 1998) , and QR = c s R/H for a thin accretion disk, we then 
have straight away 

v 2 A = a ss c 2 s . (5) 

When modeled as a mean-field — Q dynamo (c.f. Parker 1979; Reyes-Ruiz Sz Stepinski 
1995; Blackman 2000), it can be shown that due to the differential shear, the mean toroidal 
field exceeds the saturated mean poloidal field strength in the disk by a factor of H/L. 
From (|]) and (^) we have L ~ ol^H and thus 

B p ~ ajffl,. (6) 

Whether this relation holds true in the base of the disk corona where the jet would launch 
is unclear, but since the field in the corona is no greater than the field in the disk, the 
overall disk mean field provides a good upper limit at least to the total mean field energy 
density. Note that the total mean stress (B^Bp) = (B^B P ) + (b^b p ), where the second 
term is a correlation of fluctuating fields. Actually, both of these terms can drive a wind. 
Exploring this point is outside the scope of the present paper, and here we simply stick with 
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the standard approach (e.g. Blandford & Payne 1982) where the large scale fields drive the 
wind. Thus we ignore the (b^bp) term. 

The upper limit of the mean field is of order the random field, so from (|5|) we have, for 
the dominant disk field component, 

B\ ~ ATT Pd a ss c 2 s . (7) 
From the mass continuity equation of the disk we have, 

p d = M a /(Av R nHR), (8) 
where Vr is the radial infall speed. For (thin or ADAF) disks 

v R ~ a ss c s H/R ~ a ss v k (H/R) 2 , (9) 
where v\ = GM/R, the Keplerian speed. Then using (^) and (|8|) and @, we have 

o o 1/2 R 2 aL /2 (^/g)t> fc (fl)Ma 



R 2 



(10) 



This expression can be used to estimate either the magnetic field strength or Maxwell stress 
avialable in the disk for launching a wind. 



2.2. Magnetic field of the stellar shell 

The relevant stellar field is the field at the anchoring radius of the post AGB star 
from which the wind emanates. Observational interpretation suggests that the common 
envelope ejection, which would precede disk formation, removes ~ 80% of an initially 
~ 3M star (Schonberner 1993). Convective stellar models for a 3M star (Kawalar, private 
communication 2000) then tell us that outer layer will have a convective shell containing a 
mass of order 0.01M Q . We allow for a distinction between the core of the star M c and the 
mass in which the field is anchored, M#. The latter is likely to be equal to the core mass if 
the field is generated by an AGB interface dynamo (Blackman et al. 2000), but if the field 
were somehow anchored only in the thin shell, M* might be << M c . 

Since the convective shell is likely to be differentially rotating with respect to the core 
we assume that the shear rate is of order the anchoring material's angular speed, Q*. The 
shear will stretch a dynamo produced poloidal field linearly in time. The linear streching 
could operate coherently over a vertical diffusion time t d ~ Rl/r] T , where r\ T ~ Lvt is the 
turbulent diffision coefficient for turbulent velocity vt and scale L. (This would need to 
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be calculated for specific interface dynamo models, c.f. Parker 1993; Markiel & Thomas 
1999; Blackman et al 2000). Using numbers from Kawaler (private communication 2000), 
the diffusion time could be of order a few times the convective overturn time for the largest 
eddies in the post-AGB star, about 0.05yr. In any case, for linear growth of the toroidal 
field, the toroidal and poloidal fields are then related by 

~ B p *Q*t d , (11) 

which follows from the magnetic induction equation. Although the poloidal field B pif will 
likely also depend on it will also depend on other properties of the turbulence and the 
mechanism of field origin (c.f. Parker 1979, 1993; Pascoli 1997; Soker 1998; Markiel & 
Thomas 1999; Blackman et al. 2000). For an initial treatment of the problem, we simply 
consider B v * as an input condition. 



2.3. Disk inner radius 

Now that we have the stellar magnetic field, we can calculate the disk inner radius. 
Related calculations were performed by Ghosh & Lamb (1978) and Ostriker & Shu (1995). 
We estimate the inner radius by balancing the infall ram pressure from the disk with the 
magnetic pressure of the star. For the infall velocity we use the free fall velocity, as that 
corresponds to the field strength at which the disk must absolutely truncate. Thus our 
calculation represents a lower limit on the inner radius. Note for example that when the 
Alfven speed in the disk, using the stellar field strength, equals the sound speed there, 
magneto-shearing instabilities will be shut off. At this field strength the nature of the 
angular momentum transport must change (to e.g. global modes, B.Chandran, E. Ostriker, 
2000 personal communication). We assume that a suitable change occurs and thus still use 
the larger Keplerian speed for the balance below. 

Using the r~ 3 dependence for a dipole field and considering only the values at the 
equator r = R, we have 

±* M> 

where v^i is the Keplerian (~ free-fall) speed at Note that when Q*t d > 1, B"l ~ B^, 
otherwise Bl ~ B^. This is important because the angular speed couples in the former 
case but not in the latter. 

The next step is to use the mass continuity equation (||) for Ri. The result is 

Ri = BfJil + ^lrlfl"a%\HjR 4 )^Rfl\GMr X/7 M-^. (13) 
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The inner radius will be time dependent because of the time dependent accretion rate and 
angular speed of the anchoring material. We now compute the evolution of this angular 
speed. 



2.4. Angular Speed of the Star 

The disk and stellar wind luminosities depend on the angular velocity £7* of the stellar 
material where the field is anchored. As shown above, fi* is one factor which determines 
the magnetic pressure of the star, and thus R^. For small r D , decouples from 

If material is accreting onto the star from the disk at rates based on the disk formation 
models of Reyes- Ruiz and Lopez (1999), then 

/ t v 5/4 

M a ~ 6 x 10 22 g/s, (14) 

and there may be some contribution to the rotation from the accreted material. There are 
many complications in this accretion process. Here we simply assume that the material 
imparts its angular momentum to the shell (e.g. via the magnetic field) as it accretes. 
We also assume that the accreted material does not change the thickness of the star. The 
rotation speed of the relevant stellar envelope is then determined by balancing the angular 
momentum gained from the accretion with that lost from the MHD stellar wind. We have 
the following approximate equation 

M * R *^dT = -Ma^Rl + M a ^ k ,iR- - B 2 p ^rIt d . (15) 

where Qk,i is the Keplerian angular speed at the disk inner edge, and we have assumed a 
constant R*. The last term on the right hand side is due to the magnetic torque applied by 
the wind (oc Bp^B^). This equation applies when Ri corresponds to a radius at which the 
disk is spinning faster than field lines of the star at that radius. Thus it is assumed that the 
disk accretion can spin up the stellar convective envelope. There exists a regime in which 
the star spins fast enough to disrupt the disk, and angular momentum of the star is lost to 
the disk. We do not consider that regime. Note further that if Ri were greater than the 
radius at which the co-rotation speed exceeded the Keplerian speed of the star, there would 
be no accretion (Armitage & Clarke 1996). 

Re-arranging, and expanding the Keplerian term, we have 

^ + ^(MJM, + B^/M.) ^ ( GM ^Rl f \ (16) 
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In the simple limit of no disk, only the 1st and 3rd terms of fll6|) would contribute and the 
solution would be 



where we define the magnetic spin-down time scale t r 



(17) 



B^R*r D 



This quantity is the spin-down time from MHD powered rotational energy loss. For typical 
parameters, 

yr (19) 



50 



/ B pif \ 



M m / VlOkG/ V10 n cm 



10 5 sec 



As we shall see there are two regimes of interest to pursue in analytic approximation 
depending on whether t/t m is large or small. 

For the kind of disks considered here, the second term in ([16]) can be ignored 
compared with the third. Also, for the maximum accretion rates we consider 
(i.e. ~ 10~ 3 M o yr)(t/lyr)~ 5 / 4 , a mass shell of 0.01 M will at most have a 30% 
gain in mass over the lifetime, because of the form of the time dependence, so we can 
assume M* is a constant after the 1 year initial system formation period. The ODE can 
then be solved by standard integrating factor techniques, and is in fact analogous to an RL 
circuit. We find 



f2* = CExp 
where we define 



tdt' 

f in 



q(t) 



+ Exp 



* dt' 



Exp 



t dt 

f in 



q(t')dt' 



(GM c y/ 2 M a (t)Ry\t) 

RIM, 



(20) 



(21) 



Note again that we do not consider the solution before t = lyr when the accretion disk 
is still forming. For the approximate early time solution, t/lyr < t/t m < 1, we set the 
exponentials equal to 1. The result (after the 1 yr disk formation grace period) is then 



n m (t)~a,o + woU-(-&) 1/14 ), 



where 



UJo 



3 x 10 8 (GM c ) 1/2 M a0 R, 

rJmZ 



1/2 
iO 



(22) 



(23) 
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and where the subscript indicates the quantity at its initial time (t = lyr). We have 
ignored the time dependence in Ri from Q*, while including the implicit time dependence 
from M a . This is a good approximation for this regime because R l J 2 depends on V^J 7 M' 1 / 7 . 
Note that expression allows us to define the difference between initially fast (f2*o > <^o) 
and slow (f2*o < <^o) stellar rotation. 



When t/t m > 1 the exponential in (|2Cf ) begins to evolve quickly compared to the power 

1/2 

law decay of the M a R { factor inside the last integral. Factoring this out and noting that 
the product of integrals in this last term then cancel, we have 

(t) ~ q(t)t m + (K- q(t)t m ) , (24) 

where K = eQ*(t m ) + q(t m )t m (l — e). For large times (t/t m ) > 1, the first term on the right 
of (P4|) dominates and we have 

a , A (gj^|g(^ (torti , Tn>1) , (25) 

.Bp* 7 i^r D 



where the latter represents the case for which Ri depends on f2* from (|13l) . Note that the 
equations of fl25|) should be used only when it produces a value below the corresponding 
escape speed. 



2.5. Disk Wind Luminosities 

Plugging ( TUP into we have 



L dw = M dw Vl = e d ali 2 J Bp k R 2 dR ~ ^a 1 ' 2 M a (R,/ Hi)GMj R u (26) 

where indicates the outflow rate from the disk wind and V dw is the disk wind outflow 
speed. We now consider several cases which involve using the appropriate value of Ri from 
section 2.3 in (pop, then using the appropriate solution for £7* from section 2.4 in (JT3| ) for 
Ri. We will also use QTJ) for M Q . 

We first consider the case of slow stellar rotation (f^o < ^o)- The maximally luminous 
case will have the smallest disk inner radius. For a given poloidal field, the toroidal stellar 
field will be weakest for slowly rotating stars as per our discussion above, so in this case the 



disk inner radius will be the smallest. When Ri as calculated in (|T^) satisfies Ri < i?*, we 
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use Ri = R* and the disk luminosity is given by 

e d ali 2 GM c M a ( R*\ _ iQ37 e d aJ s 2 _ 2 M a o, 22 M cA / J_\ ~ 5/4 / ^/gA 



In this, and the expressions which follow, we scale a ss to 0.01, magnetic fields to 10 4 G, 
radii to 10 11 cm, diffusion times to 10 5 yr, stellar core masses to 1 M , stellar shell masses 
to 10~ 2 Mq, angular speeds to 10~ 4 Hz and mass loss rates to 10 22 g s _1 . These scalings are 
indicated by the subscripts. 

We now consider the case of an initially rapidly spinning star. The disk and star are 
strongly coupled through Ri and we must consider the different temporal regimes relative 
to the ratios t/t m and Q^td. When t < t m the first term on the right hand side of (p^) 
dominates. Then using (|T3|) for f2*To > 1, we have 



3. 8 . 9 -45 13 



e d all{GM c )iM2 fR£\^ P _ ? x 1q37 W^M^M^ f t_ \ 28 ( Rj/HA 7 erg 



For t > t m , using (|T3| ) (under the assumption that -B^* > .Bp*, or equivalently VL^Tr, > 1) 
and (|25l) in combination with (1261), we have 



3 4 



11 



-3 , 11 



— _ _ 'i?\"5" j_ 4.3 ± ( t \ 4 I RlH \ 5 ere 

~ e d a/£ {GM c )iMZBl (^J ~ 10 38 e,<_ 2 M^Mj 0i22J B| M f —J Mj^ J (29) 

In the limit that fi*T£) < 1 or B p * > B^*, from (|13|) we obtain instead, 

J_ 8-2 13 -A- £ • 2 ^45 13 

e d ail(GM c )7Mj /J^U n35 ^<-2^ 1 M a 7 0i22 / t \ » / R^JHA 7 erg 

Afo, 3 — 12 7T ~ 7 x 10 3 n — — — .(30) 

BlRj W Bl A RJ,ii \ 50jr ' V 10 / s 

For the characteristic numbers used, (EUI) does not apply because a check of Q+td from 
reveals that Q*td < 1 for t > t m for our parameter regime. 

The above formalism is valid when the magnetic luminosity is < the accretion 
luminosity GM*M C / R^. If the magnetic luminosity exceeds the disk luminosity, then the 
disk structure would be disrupted and the formalism would have to be revisited. 



2.6. Stellar wind luminosities 



The stellar wind luminosity is given by (|3|). We note again the hidden dependence of 
the field on the rotational velocity. If a dynamo generates the stellar field, then both B p 
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and B,/, can depend on For example the rotation would be required to generate the 
pseudoscalar helicity that sustains a steady field in dynamo models. Also, the toroidal 
field may induce poloidal field by springing outward. More subtle dependences should 
be considered in future work, but the important point here is the recognition of some 
dependence on the rotation, and thus on the supply of angular momentum from the disk. 
Here we simply assume that the poloidal field of the star is kept at a constant value either 
by a dynamo or by buoyancy of toroidal field, while the toroidal field depends linearly on 
the rotational speed as implied by growth from shear. In this case, equation @ gives 

L sw ~ e s B* p B^Rl ~ B^nlmRl. (31) 

We now consider the same temporal regimes for the stellar wind power that were examined 
for the disk wind. For a rapidly rotating star in the t < t m regime, the first term on the 
right hand side of (p2|) is dominant and from (|31~D we then have, 

L sw ~ t s B 2 p ^r D Rl = 10 38 e s Bl A nl 0i _,T D:5 Rl n erg/s. (32) 

For a very slowly rotating star in the t < t m regime, the second term on the right of 
( P2|) is dominant. The spin evolves to within a factor of 1/10 of uj after a few years. In this 
case f2* is such that Ri from (|13"1) is less than R* so we take Ri = R*. The stellar rotation 
rate and wind power thus approach 



O t ~ (3 x 10 7 ) 

and 



GM c Ml Q 

\ RIM* ■ 



(33) 



T Qv1 n3l B l*A TD ^ M c,l M lo,22 I , QA s 

L sw ~ 3 x 10 e s — ^ —2 — erg/s. (34) 

respectively. The slow time evolution in (|22|), which results in the approximation that the 
angular speed and luminosity are constant in (^) and (Pi|, reflects the fact that during 
this period, the angular momentum gain from accreted material and the loss of angular 
momentum from the wind nearly balance. 



For t>t m , but for VL^td > 1, using (25 ) and ( 13|) with (pl|), we have 



12 K 2 



T e s al(GM c )sMa 5 f H t y n28 e a a s 5 a ,_ 2 M £ ? 1 M a g i22 

L sw -TIK— TT = 10 13 




For t > t m , when Q*td < 1, we have instead 

L e s al(GM c )lMj /fty _ ^ esa^M^M^ ( j_ \ ^ f Rj/HA ^ erg 
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which is larger in comparison to (|35|) because a smaller field means a slower decay of Q*. 
For the characteristic numbers used, Q^td from fl25|) satisfies VL^td < 1 for t > t m , as also 
mentioned above, so here only ( PB| ) would apply for our parameter regime and not (P3|). 



2.7. Slow and Rapid Initial Stellar Rotation Cases 

In fig 1 we show plots of the disk to stellar wind luminosity ratios using e s = e<i for two 
different stellar rotation regimes. 

Slow Rot at ion: In fig la and lb plots of Ldw/ L sw for the slow rotator case (f2*o < to>o) 
are given. Fig la shows essentially the ratio of (^7]) to ( |34"D and fig lb shows the ratio (|27|) 



to (|36|). For both the t « t m and t » t m regimes, the disk wind luminosity strongly 
dominates the stellar wind luminosity by orders of magnitude. This dominance decreases 
with time in the former regime and increases in the latter regime. However the main feature 
is that for the entire range, the disk wind dominates the stellar wind. 

The star is therefore slaved to the disk throughout the evolution in this case. The 
stellar wind and spin axes are also slaved to those of the disk because the toroidal field is 
stretched perpendicularly to the axis of rotation and the wind emanates perpendicularly 
to the toroidal field. The disk and stellar winds would thus be coaxial, but the disk wind 
would strongly dominate. 

Rapid Stellar Rotation: In fig lc and Id we present plots of Ld w /L sw for the rapid 
rotator case (fi*o > ^o)- F° r t < t m the plot shows basically the ratio of (|28D to (^) and for 
t > t m the plot shows the ratio of fl30f) to (|36|) . For t < t m , the stellar wind is comparable 
to, and even slightly dominates the disk wind luminosity for the parameters used, while 
for t > t m the disk wind dominates strongly. As in the slow rotator case, the disk wind 
domination decreases slowly for t <t m and increases for t > t m . 

Thus for the large f2*o case, with f2*o near maximal speeds, the stellar wind is 
independent of the disk until t = t m . After this time its MHD wind luminosity is negligible 
compared to that of the disk. (If there were no disk at all, one would see only the rapid fall 



of of the stellar wind as in equation (|17|) ). The initial independence of the star from the 
disk means that the axes of the two winds can be misaligned. The origin of the stellar spin 
and magnetic axes are decoupled since the spin axis of the disk is independently determined 
by the initial plane of the binary orbit. Thus the disk and stellar wind system can produce 
a multi-polar wind system of up to ~ 10 38 erg/s, depending on the efficiency of conversion 
of rotational energy. 



- 15 - 



3. Discussion 

Our results show that powerful magneto-centrifugal winds can be driven in PNe 
systems. In the context of accretion disks formed from the break-up of a secondary star in a 
binary system, these winds are time dependent and have a finite lifetime. Thus our results 
would predict that transient, magnetized outflows can be driven during the transition from 
AGB to PNe. Magneto-centrifugally driven winds potentially solve two problems which 
have emerged from PNe studies (Frank 2000). First, the large powers computed above can 
explain how strong winds can be driven independently of radiation line driving processes as 
required by observations in the PPNe phase (e.g. Alcolea et al 2000). Second, they may 
explain the high degree of collimation seen in some sources (e.g. Sahai & Trauger 1998), 
and multi-polar outflows when the disk and stellar winds are contemporaneous. Below we 
discuss these points further. 

3.1. Accounting for the power of PN or PPN winds: 

The calculations above show that magnetically driven winds can account for the large 
powers observed in PPNe. Alcolea et al. (2000) make the case that observations of PPNe 
are the best place to infer the required powers and shaping mechanisms of PNe. None of 
the 3 objects for which Alcolea et al. (2000) have measured kinetic ages (OH 231.8; M 1-92; 
M2-56) can be radiation driven at the PPNe stage. The required outflow powers range from 
10 35 erg/s to 10 37 erg/s. The kinetic ages are fairly similar, ranging from 750 yr to 1800 yr. 
The radiation driving falls short of being able to provide the power by at least two orders of 
magnitude. The 9 objects for which Alcolea et. al (2000) have measured wind energies but 
not kinetic ages look reasonably similar in terms of radiation power and energetics. Since 
the kinetic ages for the few measured objects are not widely varying, it is likely that others 
in the population of 9 also probably require high outflow powers. 

The powers computed in section 2 for MHD wind driving can be as high as 10 38 erg/s 
for the parameters used. Such powers are therefore sufficient for PNe and PPNe. There is 
likely to be some efficiency factor with which the driving can take place, so to explain a 
10 37 erg/s wind, we would require an efficiency factor of e ~ 10%. More study is needed to 
determine the e factors dynamically. 

For t > t m , the disk wind power falls by two orders of magnitude or so, while the stellar 
MHD wind is of negligible power. The MHD disk wind could still supply e.g. ~ 10 36 ergs/s 
after ~ 50 yr, continuing to fall as t from (|30D . Specific examples of kinetic powers required 
in this approximate range include Hb5 (V w < 400 km/s and M ~ 6 x lO _5 M /yr; [Corradi 
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& Schwarz 1995; Pishmish et al. 2000]) thus with a kinetic luminosity L w < 3.6 x 10 36 erg/s; 
IC4406 (V w ~ 65 km/s [Corradi al. 1997] and M ~ 10~ 4 M Q /yr [Sahai et al. 1991; Cox et al. 
1991]) with L w ~ 2 x 10 35 erg/s; M2-9 (V w ~ 46 km/s [Solf 2000] and M~4x lO" 6 M /yr 
[Zweigle et al. 1997]) with L„~ 5 x 10 33 erg/s. Observations of NGC 3242, NGC 6828, 
NGC 7009, and NGC 7662 PNe (Balick et al. 1998; Dwarkadas & Balick 1998) are also 
consistent with 10 33 erg/s < L w < 10 35 erg/s. 

In all dynamo in either a disk or the star or both may be responsible for the 

field driving the wind. Magnetic flaring and the production of hard coronal X-rays might 
then be expected in PPN engines. Observations with Chandra and other X-ray telescopes 
are therefore be highly desirable. 



3.2. Accounting for the shapes of PN or PPN winds ?: 

Given that the fields can power jets, they may also play a role in shaping them. One 
possible approach to MHD collimation in these sources is simulated in Garcia-Segura et al. 
(1999). Since we have not yet dynamically calculated or simulated collimation mechanisms 
as applied to planetary nebulae the discussion that follows is speculative. The fact that the 
power of the disk and stellar winds can be comparable motivates us to proceed. 

It is known that in the context of disks (e.g. Konigl & Pudritz 2000 for review) and 
even in the context of stars (Tsinganos & Bogovalov 2000) that magnetic self collimation 
can be an important effect. The literature on collimation of disk winds is extensive. It is 
noteworthy that in emphasizing that fast rotating stellar magnetospheres can also collimate 
winds, Tsinganos & Bogovalov (2000) point out that such collimation from an underlying 
rotator will likely be effective when the quantity 



\4>eJ V fi ©/ \M C J \Vqs 

is greater than 1. Here is the magnetic flux anchored on the rotator spinning at angular 
speed Q, M w is the outflow mass loss rate, V w is the outflow speed, and these values are 
scaled to the respective values for the Sun, namely V Q = 400km/s, M Q = 1.6 x 10 12 g/s, 
S] = 3x 10~ 6 / s and ^0 = 10 22 G-cm 2 . We then have 

1/2 / N -3/2 

V c 



Q/Q Q ~ 4.2 ^ iQ26G - ^ j y-^zrj; )y 6x 10 2i g/s ) ^OOkm/sJ ' (38) 

where we have scaled to our fiducial post-AGB star parameters, taking the upper limit of the 
magnetic flux, and assuming a kinetic luminosity of the wind L kin ~ M SW V^ W ~ 10 37 erg/s. 
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We can see that Q/Q Q can be greater than 1 for the chosen parameters, and thus collimation 
in excess of that for the solar wind can occur. If we instead use representative outflow 
parameters for later stages of PNe (e.g. Balick et al. 1998), such as M = 5 x 10~ 7 M Q yr _1 
and V = 200 km s _1 ), we find Q/Qq > 100 which implies significant collimation. In general, 
these estimates suggest that self-collimation of magnetized winds in PPNe and PNe might 
be possible without requiring shocks to aid in the process. The shock-aided approach 
has been studied in Garcia-Segura et al. (1999). More study is warranted. The relative 
self-collimation of the disk wind and the stellar wind may be different for different systems. 
This needs to be studied. 

Assuming collimation can ensue, the shape of any multi-polar magnetically driven 
outflow in our picture will also depend on the stellar spin down time scale t m , and the time 
scale for the star to become a PN, = tpjy. The ratio of t m /tpN will determine how much 
shaping will be controlled by the a combined star+disk wind. When t m /tpN « 1, the 
stellar wind is likely to provide micro-structures such as knots, bullets, or ansae in mature 
PNe. When t m /tpN > 1 then structures such as multi-polar bubbles, which develop during 
the period when L dw fa L sw , may leave a large enough imprint to survive into the PN phase. 

Another important parameter is the ratio of the wind momentum densities, 
Pdw/Psw = PdwVdw/ (p sw V sw ) where the ratio to be used here is that computed for each type 
of wind independent of the actual presence of the other. The momentum density depends 
on collimation and is thus a function of angle from a wind's axis. The terminal wind 
speed of material at a given radius from any wind axis depends on the details of the wind 
launching process but should be of order several times the the escape speed at the radius 
at which the field lines threading that radius are anchored (e.g. Konigl & Pudritz 2000 
for a review). The ratio Pd w /p sw can in part determine the shape of the resulting observed 
outflow combination by determining which wind-driven "lobe" would extend farther when 
the two winds are superimposed. The difference in initial launch times of the two winds is 
also very important in determining the observed shapes. 

For two concurrent winds satisfying Pdw/Psw > 1 on average, and for which the star is 
rotating slowly {VL* < u ), we might expect two nested aligned bipolar bubbles (as shown 
in Fig 2a) when the terminal speed of the stellar wind, calculated independently from the 
disk wind, is faster over most of the outflow cross-section. The stellar wind can sweep up 
disk material, but because of the lower overall momentum density, it cannot proceed as far 
in distance. An example of a PN with this shape is Hubble 12 (Welch et al. 1999). When the 
stellar and disk magnetic and rotational axes are misaligned, a nested multi-polar bubble 
may emerge as shown in Fig 2c. A possible example of such an outflow is M2-46 (Manchado 
2000) Please note again: Fig 2 represents cartoon speculations, not simulations. 
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When Pdw/Psw < 1 and the star is rapidly rotating (Q to > u ), the stellar wind could 
push ahead of the disk wind and may produce the features seen in e.g. NGC 7007 (Fig 2b, 
Balick et al. 1998). In the case where the star is initially rapidly rotating and the stellar and 
disk axes are not aligned, then a second form of the nested multi-polar bubble may emerge 
as is shown in Fig 2d. An example of such an outflow may be He2-155 (Sahai 2000). 

In the rapidly rotating case, the stellar wind remains more or less constant over a 
spin-down period while the disk wind falls. Thus we might expect the stellar wind to 
appear more or less continuous, whereas the disk wind could produce structures only of 
comparable luminosity to the stellar wind for the first several years. In this initial phase, 
antipodal knots amidst a continuous stellar wind might appear where the disk and stellar 
winds interact, or where the disk interacts with the ambient medium. Such processes might 
have been occurring in the Egg nebula (Jura et al. 2000). 

Given that typical sound speeds for observed outflows are of order 20-30 km/s, while 
the wind velocities can be > lOOkm/s, we can estimate the size of outflow structures 
compared with the distance from the star. For a hybrid disk-star outflow with strongest 
power phase lasting ~ few years, but with an age or order hundreds of years, one should 
see knots (or shells depending on collimation) with a width of order 10 16 cm. If the outflow 
velocity were 100 km/s, then the knot/shell distance from the central core would be ~ 5 
times the shell thickness. Such features are seen in the Cat's Eye nebula (Miranda & Solf 
1992) 

Note that two interacting winds can lead to more than 2 multiple symmetry axis due to 
edge features. This and the features described above could also result from time dependent 
effects of separate interacting stellar winds or separate interacting disk winds where the 
star or disk is allowed to precess (e.g. Manchado et al. 1996). Precession would enable 
successive outbursts to be emitted on multiple axes. If there is precession, then one might 
see an additional point symmetry. These ideas may be applied to other coupled disk+star 
systems (Soker & Livio 1994, Mastrodemos & Morris 1998). 

4. Summary and Conclusions 

We have shown that MHD winds from either disks or post-AGB stars can comfortably 
power PNe outflows using reasonable estimates of stellar and disk field strengths and 
rotation parameters. The available MHD wind power turns out to be many orders of 
magnitude larger than that available for radiation/line driven outflows inferred from 
observations of PPNe (e.g Alcolea et al. 2000). Our study also highlights the importance 
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of understanding the interplay between disk winds and stellar winds in contributing to the 
diversity of observed multi-polar PNe morphologies. 

We derived maximal disk and stellar MHD wind powers for a coupled disk-star system 
at times t < t m , t > t m , (where t m is the stellar spin down time) and for fast and slow 
stellar rotators, VI^/uoq < or > 1, and values of £l*r D < or > 1. Our calculations take the 
initial time to be the onset of disk formation in a binary system from common envelope 
expulsion. The initial onset of the star and disk wind should be nearly contemporaneous. 
The layers of the star with the strong magnetic field that power the wind are only exposed 
after the common envelope is ejected, and the accretion disk forms within a year from the 
common envelope ejection. 

The resulting luminosities for both winds depend only on the disk accretion rate 
and properties of the star. This is because the inner regions of the disk are the main 
contributor to the disk wind and the disk properties can be constrained using properties 
of magneto-shearing disks and the interaction with the star. We now summarize some 
numbers for the maximum wind powers: 

* For t « t m and our choice of fiducial scalings, the maximum MHD luminosity of a 
stellar wind from an initially rapidly rotating post AGB star with field strengths consistent 
with those estimated for AGB dynamos of ~ 10 4 Gauss satisfies L sw < 10 38 erg/s, and is 
approximately constant in time. The disk wind satisfies L dw < 7 x 10 37 (t/lyr)~ L6 as the 
accretion rate falls. 

* After t m » 50yr, L sw < 10 30 (i/50yr) 2-1 , and the disk wind power goes as 
L dw < 7 x 10 35 (t/50yr)- L6 . 

* Whilst the dominance of disk vs. stellar wind depends on for t << t m , for alH >> t m , 
the disk wind always dominates, regardless of 

Since the MHD disk winds dominate the stellar winds for t > t m regardless of the 
initial stellar spin rate, the star is slaved to the disk in that regime. For t < t m , the stellar 
wind power can equal or exceed that of the disk, and in this regime the stellar wind can 
emanate along a different symmetry axis than the disk wind. Multi-polar outflows can then 
be produced. For a rapidly rotating star, the stellar MHD wind power is steady for a time 
of order t — t m , while the disk wind power falls rapidly. For t > t m , the disk wind power 
falls less rapidly than the stellar wind power. 

Generally, we would suggest that if a sustained MHD driven bipolar outflow extending 
to the core lasts for > 500 yr, the disk is dominant. The choice of scalings we use in the text 
were mainly to provide a general framework. The spin down time t m could be longer than 
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our scale of 50 yr depending on radius and stellar convective diffusion time, (c.f. eqn (|18|)). 
Disk winds imply initial formation from a binary system. An initially rapidly rotating 
star would slow down rapidly by this later time, and would only be able to contribute 
significantly if its wind were radiatively driven at this stage. If a stellar wind operates 
powerfully during t < t m it means that the central star was initially rotating rapidly. 

There is no guarantee that a binary system would form, or that the star would be 
initially rapidly rotating for every post AGB star. Thus for all systems with bipolar 
outflows, observations that can reveal the rotation rate of the central star, its magnetic field, 
or the presence of binary companions are ultimately essential for testing MHD disk/stellar 
wind paradigms. Further studies of the disk wind-stellar wind interplay are most certainly 
needed. Given the ubiquity of the accretion-central object-outflow connection throughout 
the universe the same kind of investigation may ultimately be important for other stellar 
and compact object accretion systems as well. 
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Figure Captions 

Figure 1: Plot of the MHD wind power ratio for the disk and star for four regimes with 
e s = td- a) t < t m for the initially slowly rotating star b) t > t m for the initially slowly 
spinning star c) t < t m regime for initially fast rotating star, d) t > t m regime for initially 
fast rotating star. See text. 

Figure 2: Cartoon representing our speculation on a possible correspondence between 
between HST PN images and multi-polar wind structures from the hybrid stellar-disk wind 
paradigm, for different initial stellar spins and thrust ratios Ud w /n sw (see text). (Note that 
these are not simulations, but speculations.) a) top left: < uoq and Tldw/n sw > 1 so star 
is slaved to disk and the only possibility is nested winds. Inset is PN Hubble 12 (Welch et 
a/,1999). b) top right: Q* < and U dw /U sw < 1. Inset is NGC 7007 (Balick et a/.1998). 
c) bottom left: f2*o > and U dw /Ils W > 1- Inset is M2-46 (Manchado et a/.2000). d) 
bottom right: f2* > and Tld w /Tls W < 1- Inset is He2-155 (Sahai 2000). See text. 
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